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ABSTRACT 



We discuss interstellar temperature determinations using the excitation equi- 
librium of the 2 P levels of Si II and C II. We show how observations of the 2 -P3/2 
fine structure levels of Si II and C II (which have significantly different excitation 
energies, corresponding to ~ 413 and 92 K, respectively) can be used to limit gas 
kinetic temperatures. We apply this method to the z = 4.224 damped Lyman-a 
system toward the quasar PSS 1443+27. The lack of significant absorption out of 
the Si II 2 P3/2 level and the presence of very strong C II 2 P3/2 provides an upper 
limit to the temperature of the C II*-bearing gas in this system. Assuming a 
solar Si/C ratio, the observations imply a 2a limit T < 954 K for this absorber; a 
super-solar Si/C ratio gives stricter limits, T < 524 K. The observations suggest 
the presence of a cold neutral medium; such cold gas may serve as the fuel for 
star formation in this young galaxy. 

Subject headings: ISM: evolution - ISM: atoms - galaxies: evolution - quasars: 
absorption lines 



1. Introduction 

High-redshift damped Lyman-a systems (DLAs) are the highest column density class of 
QSO absorption lines. Defined by log A(H I) > 20.3 (Wolfe et al. 1986), these systems are 
thought to trace the interstellar medium (ISM) of high-redshift galaxies. Dedicated surveys 
over the past two decades have helped trace the global properties of high-redshift DLAs, 
including their contribution to the cosmological baryon density (Storrie-Lombardi & Wolfe 
2000; Prochaska & Herbert-Fort 2004), their chemical enrichment (e.g., Prochaska et al. 
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2003), their dust content (e.g., Pettini et al. 1994), and molecular fraction (Ledoux, Petitjean, 
& Srianand 2003). These studies have demonstrated that the DLAs have a baryonic mass 
density comparable to the mass density of modern galaxy disks, that the metallicity of high- ,2 
DLAs is slowly increasing, and that the majority of DLA sight lines have low dust-to-gas 
ratios and molecular fractions. 

Understanding the detailed physics of the ISM in DLAs is an important step in un- 
derstanding high-redshift galaxies in general. Wolfe, Prochaska, & Gawiser (2003a) have 
constructed detailed models for the thermal equilibrium of the ISM in a set of DLAs, calcu- 
lating the heating rate experienced by the gas due to the ultraviolet emission from young hot 
stars. They used observations of absorption out of the 2 P^/2 level of C II (hereafter C II*) - a 
direct indicator of the cooling rate through [C II] 158 /xm emission (Pottasch, Wesselius, van 
Duinen 1979) - to infer the actual heating rate experienced by the gas (assuming thermal 
equilibrium). Their comparison of the inferred and calculated heating rates suggests DLAs 
harbor significant star formation. 

The models of Wolfe et al. (2003a,b) required that the observed C II* in DLAs arise 
in a cold neutral medium (CNM), i.e., in gas with temperatures T < 1000 K; their WNM 
models give SFRs that violate observations of the bolometric background. The detection of 
H2 absorption in some high-z DLAs is further evidence that at least some of these systems 
contain a CNM (Ledoux et al. 2003; Hirashita & Ferrara 2005). However, such temperatures 
are at odds with 21-cm absorption studies; all z > 3 DLAs searched for 21-cm absorption 
show T s > 1400 K (2(j) (Kanekar & Chengalur 2003). Wolfe et al. (2003b) argue that this 
discrepancy is likely due to the differing properties of the sight lines probed by the optical 
background sources and by the more extended radio sources. 

In this Letter we present a method for determining the kinetic temperature of interstellar 
matter based solely on basic atomic physics. Our method compares the excitation of the 
upper 2 Pz/2 fine-structure levels in Si II and C II, which have excitation energies that differ 
by a factor of four. We describe our approach in §2. We apply this technique to the 
z abs ~ 4.224 DLA toward the quasar PSS 1443+27 in §3, demonstrating that this DLA 
contains a substantial reservoir of cold gas. Lastly, we discuss the implications of this 
temperature determination in §4. 

2. Determining Gas Temperatures from iV(Si II*) and N(C II*) 

In this section we discuss how analysis of the 2 P fine structure excitation in Si II and 
C II can be used to limit interstellar kinetic temperatures. Srianand & Petitjean (2000) have 
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used a similar analysis to limit temperatures and densities in an associated absorber. Silva 
& Viegas (2002, hereafter SV) give a detailed summary of the excitation equilibrium of fine 
structure lines in DLAs. 

We treat the excitation of the 2 P fine structure states as a two- level atom. 1 The equilib- 
rium ratio of the densities in the upper and lower levels (Spitzer 1978), n 2 and ni, respectively, 
is 

B 12 U„ 12 (z) + F 12 + Yl n kll2 

^1 = * . (1) 

ni A 21 + B 2 iU Ul2 (z) + T 21 + n k l 21 

k 

The quantities A 21 , B 21 , and B u are the familiar Einstein transition probabilities. The 
energy density of the cosmic microwave background, u Ul2 (z), for direct excitation of the 
transitions is calculated assuming a standard cosmology, i.e., Tqmb = Tq(1 + z), where 
T = 2.725 K (Mather et al. 1999). The fluorescent rates T 12 = T 21 = due to the great 
opacity of the ground-state transitions of interest (Sarazin, Rybicki, & Flannery 1979; see 
Wolfe et al. 2003a). The summations describe excitations and deexcitations with collision 
partners k(= e~ , p + , and H°), where n k is the particle density of each partner and 7^ and 
7^ are the Maxwellian-averaged collision rate coefficients for excitation and deexcitation. 
The collision rates are related to one another and to f2^ 2 (T), the collision strength (Spitzer 
1978), by: 

7 f 2 oc 7* exp (-kT 12 /kT) oc Q k 12 (T) T~ 1/2 exp (-kT 12 /kT), (2) 

where kT 12 is the energy of the 2 Ps/ 2 level above the ground state, and T is the gas kinetic 
temperature. We adopt atomic data from the same sources as SV throughout. 

The column density of material in the 2 Pz/2 level of C II is 

TV(Cir) w f^^A(C)n H ds, (3) 
J n c 11 

where A(C) is the gas-phase abundance of carbon, hh is the density of neutral hydrogen, ds 
is the differential pathlength along the line of sight, and the ratio ^cn*/ n cn is a function 
of four quantities: z a b s , nn, x = n e /ri}i ~ n p /nn, and T (Eqn. 1). We have assumed 
riQYi n c = A(C)n H , consistent with limits on C I in DLAs (see Figure 12 of Wolfe et al. 
2003b). A similar expression applies for Si II*. The column densities of C II* and Si II* arc 
therefore density- and excitation-weighted integrals over pathlength. 

To understand the usefulness of Si II* and C II* to limit gas temperatures, imagine that 
the excitation of the 2 Pz/ 2 states were only due to electron collisions and deexcitation due to 
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spontaneous emission. In this case iV(Si U*)/N(C II*) oc [A(Si)/A(C)] x [n e 12 (T) Si /n e 12 (T) c ] x 
exp{ — (T^ 1 — T^)/T}, i.e., the ratio is only a function of temperature and the Si/C ratio. 
Figure 1 shows the ratio N(Si U*)/N(C II*) as a function of temperature, including all of the 
terms in Eqn. (1), for several densities at z ~ 4.2 assuming [Si/C]=0. 2 iV(Si II*)/A r (C II*) 
varies slowly with temperature for T > 1, 000 K (for all but the highest densities) and more 
strongly at lower temperatures. The sensitivity to density is modest, especially for low x, 
because the populations of the 2 P 3 / 2 levels of Si II and C II depend on density in a similar 
manner. We use the ratio of the 2 P%/2 to 2 Pi/2 levels of C II and Si II to constrain the 
allowable range of densities. 

Gas at T < 1,000 K produces ratios logiV(Si U*)/N(C II*) < -2.8; the corresponding 
optical depth ratio of the 1264.738 A line of Si II* to the typically-unresolved 1335.663 and 
1335.708 A doublet of C II* is r(1264)/r(1335) < 0.015 (/-values from Morton 2003). Thus, 
the use of this ratio to distinguish cool from warm temperatures requires large C II* optical 
depths and good S/N at Si II*. Observations of iV(Si U*)/N(C II*) will allow a range of 
parameters nn,x, and T, which apply to the C II*-bearing gas. However, the maximum 
allowable temperature can be low if the limits on JV(Si II*)/jV(C II*) are low. 

3. Application to the z = 4.224 Damped Lyman-a System Toward 

PSS 1443+27 

Prochaska et al. (2001) present high- resolution (R ~ 40, 000), high S/N observations of 
the z = 4.224 DLA toward of the quasar PSS 1443+27 using HIRES (Vogt et al. 1994) on 
the Keck I 10-m telescope. The properties of this DLA are summarized in Table 1. Due to 
its redshift and relatively strong metal lines, column densities can only be measured for a few 
species. Table 1 gives a measurement of Fe II (from the 1611.2 A transition) and meaningful 
limits on C II*, Si II and Si II*. All but the last of these are from Prochaska et al. (2001). 

Figure 2 shows the HIRES spectrum at the expected location of Si II* in this absorber, 
including our fitted continuum. Our limit to Si II* was determined using an empirically- 
measured signal-to-noise ratio (following Sembach & Savage 1992) and includes continuum 
placement uncertainties. Figure 3 shows the normalized profiles of Si II* and C II*. Also 
shown are models for the expected Si II* absorption for a canonical warm neutral medium 
(WNM; with n H = 1 cm' 3 , x = 0.1, and T = 8,000 K) and CNM (with n H = 10 cm" 3 , 
x = 10~ 3 , and T = 80 K). The models are derived by scaling the C II* profile by the 
amount predicted using Eqns. (1) and (3) and assuming [Si/C]= 0. The hypothesis that 
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C II* absorption toward PSS 1443+27 arises in a canonical WNM is not consistent with the 
observations. 

To limit T, we calculate the excitation balance of the fine structure levels of Si II and 
C II for densities —4 < logn e < +2 and —2 < logn H < +3 and temperatures 1.5 < logT < 
4.5. We compare these calculations with our observational limits for N(Si U*)/N(C II*), 
iV(C H*)/JV(C II), and iV(Si H*)/JV(Si II) to constrain the physical properties (n H , x, and 
T) of the gas in the C II*-bearing gas. 3 

The allowable range of physical conditions in the PSS 1443+27 DLA then depends upon 
the ratios of Si/C and C/Fe, which are used for determining the expected iV(Si U*)/N(C IP) 
and N(C H*)/N(C II), respectively. C II is rarely measurable in DLAs due to satura- 
tion effects; we estimate N(C II) from iV(Fe II) (Table 1) with an assumed abundance 
ratio. As a fiducial, we adopt the solar system ratios from Grevesse & Sauval (1998): 
log(C/Fe) = +1.02 and log(Si/C) = -0.96. This gives logiV(C II) = 16.22 + 0.06. Nucle- 
osynthetic effects are unlikely to modify the C/Fe ratio, as a roughly solar value is found in 
low-met allicity Milky Way stars (e.g., Carretta, Gratton, & Sneden 2000). However, nucle- 
osynthetic effects could enhance a element abundances, leading to high Si/C ratios (see, e.g., 
Prochaska 2004). In one case where the C II measurements may be reliable, D'Odorico & Mo- 
laro (2004) find logiV(Si U)/N(C II) = -0.74+0.07 and logiV(C II)/JV(Fe II) = +1.25 + 0.12 
(assuming the Fe II column from Lu et al. 1996). This implies [C/Fe] = +0.23 and 
[Si/C]= +0.22 compared with our adopted solar system abundances. The overall metallicity 
of this absorber is [Si/H] = —1.72, however, which is significantly lower than the metallicity 
of the DLA studied in this work. 

The presence of dust will raise the C/Fe ratio due to differential depletion. The deple- 
tion of Fe in DLAs is typically modest, similar to the "halo" and "disk-halo" clouds in the 
Milky Way (Savage & Sembach 1996) with [Fe/H]> -0.9. 4 For diffuse clouds in the Milky 
Way, Sofia et al. (2004) find [C/H]~ -0.3. Thus, we expect [C/Fe]< +0.6. In low-depletion 
Milky Way clouds, and by extension DLAs, Si is likely not depleted much differently from C, 
and we assume [Si/C]= 0. The small number of species for which measurements are possible 
in the DLA toward PSS 1443+27 make it difficult to judge the degree of Fe depletion. Mod- 
erate depletion is suggested by [Si/Fe] > +0.17, but this could also result from a-enhanced 
abundances. 



3 These three physical quantities are coupled in the calculations. A wide range of T-dcpendcnt x and «h 
values are allowable, and only T is stringently constrained in our analysis. 
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We give temperature limits for the C IP-bearing gas in the z ~ 4.224 DLA toward 
PSS 1443+27 for six different abundance assumptions in Table 2. The first three assume 
[Si/C] = with [C/Fe] = +0.0, +0.3, and +0.6, allowing for modest levels of depletion. The 
next three assume [Si/C]= +0.3, i.e., a-enhanced abundances (see Prochaska 2004). We 
note that the model in which [Si/C] = [C/Fe] = is not self-consistent: scaling both Si II and 
Fe II by solar abundances gives inconsistent values for the C II column density. 

The highest limit is T < 954 K, assuming assuming [C/Fe] = +0.6 and [Si/C] = 0. 
This is likely an extreme assumption for the C/Fe ratio. The highest temperature limit 
for a-enhanced abundances, which we prefer due to the distribution of [Si/Fe] at very low 
metallicities (where dust depletion effects are likely small; see Prochaska 2004), is T < 524 K. 
Temperatures above 1400 K, consistent with T s limits for the z > 2.9 DLAs observed at 21-cm 
(Kanekar & Chengalur 2003), require [C/Fe] > +0.75 to +0.8 for both solar and a-enhanced 
abundances. We feel this is inconsistent with our knowledge of differential depletion in the 
Milky Way and DLAs. 

4. Discussion 

We have presented a method for limiting gas temperatures in the ISM of galaxies through 
measurements and analysis of the 2 P3/2 fine structure levels of Si II and C II. We have applied 
this method to limit the properties of the C IP-bearing gas in the z = 4.224 DLA toward 
PSS 1443+27. We rule out the hypothesis that the C IP absorption arises in a WNM. Our 
conservative temperature limit for this gas is T < 954 K; we obtain stricter temperature 
limits if [Si/C] > or [C/Fe] « 0: T < 524 K. The detection of a CNM (Wolfire et al. 1995) 
in a high-redshift DLA is significant: while the gas seen in this DLA is likely not associated 
with the dense star-forming clouds, our result demonstrates that the physical conditions of 
the ISM in this system do not preclude the existence of cold material, including, in principle, 
the dense clouds from which stars could form. 

The detection of a CNM in the z=4.224 absorber toward PSS 1443+27 is the first 
detection at such a large redshift. We stress that our measurements allow a WNM as part of 
a multiphase medium toward PSS 1443+27, but the majority of the C IP cannot come from 
warm material. 5 The existence of CNM material may be a feature of many high- z DLAs, as 
suggested by Wolfe et al. (2003a,b). While z < 2 H I 21-cm absorption-line measurements 
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suggest the presence of cold H I in DLAs, no 21-cm absorption has been found in DLAs with 
z > 3 (Kanekar & Chengalur). 6 We note, however, that at least three z > 2.5 DLAs show 
H 2 absorption (Ledoux et al. 2003), indicating the presence of cold gas (Hirashita & Ferrara 
2005). Searching for H 2 toward PSS 1443+27 might be difficult due to the strength of the 
Lyman-a forest at z ~ 4. In the future we will apply our method for determining gas kinetic 
temperatures to DLAs for which 21-cm and H2 measurements have been attempted. 

We note that the DLA toward PSS 1443+27 may not be typical. It shows a high 
metallicity ([Si/H]> —1), especially compared with other z > 4 DLAs (Prochaska et al. 
2003). Furthermore, we chose this DLA for this experiment because it has very strong C II*. 
It has the highest C IP/H I ratio measured, some +0.6 dex above the next highest (the 
z = 1.92 DLA toward Q2206-19; Wolfe et al. 2004). The intensity of radiation in this DLA 
calculated following Wolfe et al. (2004) is higher than all other systems. Although the fuel 
for star formation may be present in this DLA, no optical counterpart to this DLA has 
been identified in deep ground-based and Hubble Space Telescope images of PSS 1443+27 
(Prochaska et al. 2002; L. Storrie-Lombardi, private communication). 

The method presented here will be discussed further, and temperature limits given for 
a larger number of DLAs, in Howk et al. (2005). We note that the use of the relative popu- 
lations of the 2 P3/2 levels of Si II and C II may also be useful for constraining temperatures 
in Milky Way gas. 
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Table 1. PSS 1443+27 Damped Lyman-a System Properties 51 



log AT (Si II*)/iV(C II*) < -3.47 (2a) 



a iV(H I), [Fe/H] and column densi- 
ties, with the exception of Si II*, from 
Prochaska et al. 2001. 



Table 2. Temperature Limits for z = 4.22 Damped System 



Model 


[Si/C] 


[C/Fe] 


T[K] 


logn H a 


1 


+0.0 


+0.0 


< 478 


1.44 


2 


+0.0 


+0.3 


< 724 


1.04 


3 


+0.0 


+0.6 


< 954 


0.56 


4 


+0.3 


+0.0 


< 416 


1.46 


5 


+0.3 


+0.3 


< 436 


1.14 


6 


+0.3 


+0.6 


< 524 


0.64 



a Hydrogen density giving maximum al- 
lowed temperature in each model. 



Quantity 



Value 



logiV(Si II) 
logiV(Fe II) 
logiV(C I) 
logiV(C II*) 
logiV(Si II*) 



log7V(H I) 
[Fe/H] 



4.2240 
20.80 + 0.10 
-1.10 + 0.11 

> 15.43 
15.20 + 0.06 
13.37 + 0.09 

> 14.71 
< 11.24 (2a) 



2.0 2.5 3.0 3.5 4.0 

log T 



Fig. 1.— Predicted JV(Si E*)/N(G II*) for a DLA at z ~ 4.22 with logn e = -1 and 
a large range of hydrogen densities. The dashed line represents the 2a upper limits to 
N(Siir)/N(CU*) for the DLA toward PSS 1443+27. Electrons generally dominate the 
collisional excitation of Si II and C II (7^ S> 712)- Although at very low values of x = rie/^H 
(bottom curves), H I collisions become important. The change in slope of the lowest-x curves 
is due to the larger relative importance of H I in collisional excitations of C II compared 
with Si II (compare Figures 3 and 7 of SV) and the differing temperature dependence for 
7^. The ratios of the 2 P3/2 to 2 Pi/2 levels of C II and Si II can be used to constrain n H and 
x. 
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Fig. 2. — Observed spectrum near the expected position of Si II* for the z = 4.224 DLA 
toward PSS 1443+27 (shown by the dashed line). The thick grey line represents our fitted 
continuum. The dotted region denotes the location of two poorly-subtracted sky emission 
lines, which cause the narrow depressions in this region of the spectrum. 



13 



1.10 F 




-100 -50 50 100 
Velocity (km/s) 

Fig. 3.— Normalized Si II* and C II* profiles for the z = 4.224 DLA toward PSS 1443+27. 
Included in the panel showing Si II* are two predicted Si II* profiles derived from the observed 
C II* profile assuming standard WNM (blue, with tir = 1 cm -3 , x = n e /n^ = 0.1, and 
T = 8,000 K) and CNM (red; with n H = 10 cm" 3 , x = 10~ 3 , and T = 80 K) properties. 
These models assume [Si/C] = 0. The arrows denote regions of definite saturation in the C II* 
profile. The dotted region in the Si II* profile denotes the location of two poorly-subtracted 
sky emission lines. 
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